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Abstract. Elliptic iV-soliton-type solutions, i.e. solutions emerging from the application of 
N consecutive Backlund transformations to an elliptic seed solution, are constructed for all 
equations in the ABS list of quadrilateral lattice equations, except for the case of the Q4 
equation which is treated elsewhere. The main construction, which is based on an elliptic Cauchy 
matrix, is performed for the equation Q3, and by coalescence on certain auxiliary parameters, 
the corresponding solutions of the remaining equations in the list are obtained. Furthermore, 
the underlying linear structure of the equations is exhibited, leading, in particular, to a novel 
Lax representation of the Q3 equation. 



1. Introduction 

In a series of recent papers, soliton type solutions of two-dimensional integrable quadrilateral 
lattice equations were explored, [6J [TOj [23]. Integrability here is understood in the sense of 
the multidimensional consistency property, |32[ I13j . which states that such equations can be 
consistently embedded in a multidimensional lattice, and which has been accepted as a key 
integrability characteristic. Several examples of equations exhibiting this property were known 
for more than two decades, cf. e.g. [361 EH US EQj 12I]> but it was only recently that a full list 
of scalar multidimensionally consistent equations was obtained, cf. [2J, and also [3]. Apart from 
lattice equations of KdV type which had been established early on, a number of novel equations 
arose from this classification, for which no further structures (such as Lax pairs, inverse scattering 
scheme or direct linearization treatments) were a priori known. As was demonstrated in [131 [23] 
some of these structures, in fact, follow from the multidimensional consistency of the equation 
itself. However, the construction of explicit solutions remained an open problem, and this 
was systematically undertaken in the series of papers mentioned above, and they reveal some 
surprising new features. In [21] we revealed some of the underlying structures of the equations 
in the ABS list, showing that they are deeply interrelated not only through degeneration, but 
through Miura type relations as well. In particular, closed-form N solitons for the whole ABS list 
were obtained, with the exception of the "top" equation in the list, the so-called Q4 equation, 
which was first discovered by V. Adler, [JJ, as the permutability condition for the Backlund 
transformations of the famous Krichever-Novikov equation. Recently we have given explicit 
7V-soliton solutions for Adler's equation [11] using a new constructive approach in which the 
solutions emerge in Hirota form. With respect to the Cauchy-matrix approach considered here, 
the equation Q4 has special features which warrant a separate treatment. The main focus in 
[24] . cf. also [7], was on the construction of "rational soliton solutions" for the Q3 equation, 
which turns out to be the equation in which the features of the entire remaining list culminate. 
Its underlying structure incorporates ingredients from all the other equations in the ABS list, 
and consequently from its explicit iV-soliton solutions, the corresponding solutions of all the 
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"lower" equations in the ABS list are obtained through degenerations (i.e. via limits on some 
fixed parameters). Thus, we were able to present explicitly all such rational N-soliton solutions 
for those equations in closed form, (cf. also [19] for bilinear form and Casorati determinant 
expressions of the same solutions). 

In this paper we will generalise the Cauchy-matrix approach of [24] to the elliptic case, i.e. 
to solutions which we could call elliptic N-soliton solutions. The problem of finding such elliptic 
solutions should not be confused with the one of finding soliton solutions of the elliptic equation 
Q4 (in which the lattice parameters take values as points on an elliptic curve) considered in 
[llj . The treatment here follows closely the derivations in [24], but the lifting of the rational 
solutions to the elliptic case introduces some new interesting features, such as the emergence of 
a non-autonomous Cauchy kernel that incorporates the core strucuture behind these solutions. 
Furthermore, some aspects become actually more transparant in the elliptic soliton the 
elliptic addition formulae at the heart of the development in some sense drive the construction. 
Moreover, we unravel in this paper some of the underlying linear structures behind these solu- 
tions, which as a byproduct leads to the construction of a novel Lax pair for Q3, which seems 
more natural than the Lax pairs that are obtained from multidimensional consistency following 
the recipe in [23\ I13j . Thus, we expect that this new Lax representation can prove useful in 
the construction of wider classes of solutions, such as finite-gap solutions, and inverse scatter- 
ing solutions with radiation. We also envisage that the structures revealed in this paper may 
prove important in understanding the general nature of elliptic solutions of integrable discrete 
equations and the corresponding many-body systems of Calogero-Moser and Ruijsenaars type, 
cf. also [271 [22]. 

2. Preliminaries: the ABS list 

First, we need to establish some notations that we will employ throughout the paper. Equa- 
tions within the class of quadrilateral partial difference equations (PAEs) have the following 
canonical form: 

Q Pi q (u,u,u,u) = , 

where we adopt the short-hand notation of vertices along an elementary plaquette on a rectan- 
gular lattice: 

Schematically these are indicated in Figure [TJ Here p, q are parameters of the equation, related 
to the lattice spacing, so each parameter is identified with a direction in the lattice. It is 
useful, and actually crucial for the classification problem as solved in [2], to consider the lattice 
parameters p, q as points on an algebraic curve, i.e. given by p = (p, P) , q = (q, Q) , where the 
coordinates p, P and q, Q respectively are related through some polynomial equation. 

ABS, in [2], considered the case where u is a scalar (i.e. single-field) quantity, and restricted 
themselves to the affine-linear case, i.e. the case in which the function Q is assumed to be affine 
linear in each of its four arguments. They required the multidimensional consistency property, cf. 
[321113] . That is, in each pair of directions of the multidimensional lattice a copy of the equation 
(chosing values of the lattice parameters p, q, etc., associated with those directions) can be 
imposed such that the iteration of well-chosen initial values leads to a unique determination of 
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Figure 1. Arrangement of the shifted dependent variable on the vertices of a 
quadrilateral and association of the lattice parameters to the edges. 



the solution in each lattice point, avoiding possible multivaluedness that would generically occur 
in the iteration of the solution by evaluating the values of u on lattice points following different 
configurations of quadrialterals. In the integrable case, such multivaluedness does not occur by 
the precise structure of the equation and the combinatorics of the evaluation process. More 
concretely, this means that any three copies of the equation with different parameters, namely 



Q p ^(u,u,u, u) 



, Qp ir (u, U, U,u) = , Qq, r (u, U, U, u) = 

can be simultaneously solved, and lead to a unique and single- valued determination of the triple- 
shifted point u at the vertex of an elementary cube on which initial conditions are given by u, 
u, u, u, where the — denotes the elementary shift in an additional lattice direction associated 
with the lattice parameter r. 

The classification in [2] yielded, among others, the following list of equations (presented in a 
slightly different notation) 



Q-list: 



Ql 
Q2 

Q3 
Q4 



p(u — u)(u — u) — q 2 (u — u)(u — u) = 5pq (p — q) 

p(u — u)(u — u) — q(u — u)(u — u) + pq(p — q)(u + u + u + u) 

= pq{p-q){p -pq + q ) 
p(l — q 2 )(uu + uu) — g(l — p 2 )(uu + uu) = 



: y 2 



°2n 



) I (uu + uu) + 5 



.(1 £2) 



Apq 



p(uu + uu) — q(uu + uu) 
pQ-qP 



1 



on oo 

p £ q L 



(uu + uu) — pq(l + uuuu) 



(2.1a) 
(2.1b) 

(2.1c) 
(2.1d) 



with 5 being a fixed parameter, and where the form of Q4 as given in (|2.1d|) was actually due to 

o o 

[18] . and involves additional parameters P, Q related to p, q respectively via the Jacobi elliptic 
curve relations: P 2 = p 4 — jp 2 + 1 , Q 2 = q 4 — ^q 2 + 1 , in which 7 denotes the modulus of this 
curve. 
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In what follows we will focus mostly on Q3 and its degenerations, and for the purpose of our 
treatment we will rewrite (I2.1cp in a different parametrisation which is more adapted to the 
actual structure of the solutions, namely as follows: 

Q3 : P(uu + uu) — Q(uu + uu) = (p 2 — q 2 ) (uu + uu] + — — (2-2) 



) |juu + uuj + pg 

in which, following [23], the lattice parameters p, q are related to the parameters p, q by the 
relations 

P = -j n , q = -s 2 ' ( 2 - 3 ) 

p z — a z q 0- 

and where (p, P) and (q, Q) are points on a Jacobi elliptic curve (different from the one mentioned 
above associated with the Q4 equation) with branch points ±a, ±6, i.e. P 2 = (p 2 — a 2 )(p 2 — b 2 ) , 
and Q 2 = (q 2 — a 2 )(q 2 — b 2 ) . The fixed parameter A will have a special significance in terms of 
certain arbitrary ceofficients in the solutions, as we shall see. Even though Q3, as given in (|2. lc|) . 
has a rational parametrisation, it turns out that for the solution structure, both here (where we 
are dealing with elliptic solutions) as in the case of rational solitons as in [23], this new Jacobi 
elliptic curve plays a crucial role. Thus, throughout this paper we will consider the parameters 
p and q, and the associated parameters P, Q on the elliptic curve, to be the natural parameters 
for the solutions, and consequently we prefer to express the solutions of the degenerate cases, 
such as Q2 and Ql in terms of these parameters. Other degenerate cases involve the so-called 
H-equations from [2j, and they are given by 



H-list: 



HI: (w - w)(w - w) = p 2 - q 2 (2.4a) 
H2 : (w - w){w - w) = {p 2 - q 2 ){w + w + w + w) - p 4 + q A (2.4b) 

- - P 2 - Q 2 

H3 : P{ww + ww) - Q(ww + ww) = 2A — — — (2.4c) 

• {J 

where we have now expressed them in terms of the new parameters p and q, whereas in (|2.4c|) we 
have introduced associated parameters P, Q defined by the relations P 2 +p 2 = a 2 , Q 2 +q 2 = a 2 . 
We have omitted from the list the equations denoted by Al and A2 in [2], which are equivalent to 
Ql and Q3 respectively upon gauge transformations. In a more recent paper, [3], ABS achieved 
a somewhat stronger classification result, still in the scalar and affine-linear case, but otherwise 
under less stringent assumptions. The various equations are connected to each other through 
degenerations of the parameter curve, from Q4 down to all other equations. Concentrating 
on the Q-list alone, the corresponding coalescence diagram is given in Figure [2j The question 
as to what is known about explicit solutions of these equations is an interesting one. In fact, 
most results so far are obtained for the equations of KdV type, and go already back several 
years. Thus, identifying this subclass of equations, we note that HI is the so-called lattice 
potential KdV equation, which also appeared under the guise of the permutability condition 
of the Backhand transformations of the (continuous) KdV equation, [36], or as the so-called e- 
algorithm in numerical analysis, due to Wynn, cf. [38]. The equation (Ql)o, i-e. the Ql equation 
with fixed parameter 5 = 0, is the lattice Schwarzian KdV equation, or cross-ratio equation, 
first presented as an integrable lattice equation in [25]. The equation (H3)o can, up to a point 
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Figure 2. Coalescence diagram for equations in the ABS Q list (this was first 
given by Adler and Suris in [1]). 



transformation, be identified with the lattice potential modified KdV (MKdV) equation, also 
referred to as the Hirota equation, and (Q3)o can be obtained from an interpolating equation that 
was first presented in [30] . and which hereafter we refer to as the NQC equation following [35] . 
For all these equations multi-soliton solutions are implicit in the direct linearisation approach 
developed for lattice equations in [30\ I34| . In fact, inverse scattering type solutions including 
radiation are in principle (i.e. on a formal level) included in that scheme as well. Finite-gap 
solutions of those equations are the ones that are connected to the periodic reductions of those 
lattice equations, i.e. the reduction to finite-dimensional integrable (in the sense of Liouville) 
dynamical mappings, cf. [33]. Those were developed for the lattice KdV in [26] . whereas finite- 
gap solutions for the Hirota equation and the lattice Schwarzian KdV were constructed from a 
geometric perspective in [12|. 117]. Finally, scaling symmetric solutions, leading to reductions to 
discrete Painleve equations were constructed in [29^ 131] . To our knowledge this exhausts the 
list of more or less explicit solutions to the lattice equations of KdV type. For any of the new 
equations in the ABS list, such as the deformed equations depending on the parameter 5, as 
well as H2, Q2 and Q4, no explicit solutions were presented until more recently. The case of Q4 
was considered in [6] , whilst with the exception of Q4, all of the other new ABS equations were 
shown to be Backlund related to the older equations [8], which percipitated the generalisation of 
the Cauchy-matrix approach to these new systems and hence the results of [7] 1241 [19] . The aim 
of the present paper is to add a new class of solutions, namely elliptic soliton type solutions, of 
the equations in the ABS list to those found in recent years. 

The outline of the remainder of the paper is as follows. In section 3 we show that additiona 
formulae for elliptic functions (in the Weierstrass class) have a natural interpretation in terms 
of partial difference equations. Thus, elliptic functions are shown to arise naturally as seed 
solutions of a number of such equations. In section 4 we will set out the ground-work for our 
construction of elliptic iV-soliton solutions, defining the basic objects and relations in terms of 
elliptic Cauchy matrices. These results will subsequently be used in section 5 to arrive at the 
main statement (Theorem 5.1) presenting the elliptic iV-soliton solution for Q3. Furthermore, we 
present a realisation of those solutions in terms of a novel Lax pair for Q3, identifying explicitly 
the eigenfunctions. In section 6, we demonstrate that the closed form solution of Q3 of the 
previous section arises also from a Backlund chain. Finally, in section 7 we show how the 
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solutions of Q3 degenerate to analogous solutions of the lower members of the ABS list. Some 
direct verifications of the resulting solutions of the Q- and H-equations are also presented. 

3. Elliptic addition formulae as solutions of PAEs 

In this section we will show that many of the equations in the ABS list admit elementary 
solutions in terms of elliptic functions. In fact, one can consider the lattice equations as models 
for basic addition formulae of elliptic functions. We will then, in subsequent sections, generalise 
these elementary solutions to more general elliptic solutions, which we will call elliptic solitons, 
and we will show that the elementary solutions play the role of the seed solutions for those 
elliptic solitons as viewed as arising from a Backlund chain. 

3.1. Elliptic addition formulae. In this paper we will work throughout exclusively in terms of 
Weierstrass functions, and we employ primarily the basic addition formulae for these functions. 
In what follows, a(x) = a(x\2uj,2uj'), ((x) = C(x\2uj,2uj'), = a(x\2u, 2u/), denote the 

standard Weierstrass a-, £- and p-functions with simple periods 2cj, 2u/ £ C, see e.g. [5]. Let 
us now recall the following elliptic identities: 

(j-function: The celebrated three-term identity for Weierstrass cr-function 

a(x + a)a(x — a)a(y + b)a(y — b) — a(x + b)a(x — b)a(y + a)a(y — a) 
= a[x + y)cr(x — y)cr(a + b)a(a — b) (3.1 
can be rewritten in terms of 

a(x + k) 



<f> K {x) :-- 

a{K) a(x) 

as: 

<5>K(x)<S>\(y) = <S> K+ \(x)<P x (y -x) + $ K (x - y)®K+\(y) ■ 
("-function: The Weierstrass (-function C( x ) = cr'(x)/a(x) obeys 

a{x + y) a(y + z) a(x + z) 



((x) + ((y) + ((z) -((x + y + z) 



a{x) a(y) a(z) a(x + y + z) 
or: 

$ K (x)$ K (y) = <$> K (x + y) [C(«) + ((x) + ((y) -((k + x + y)} 
p-function: The Weierstrass p-function p(x) = —£'(x) obeys: 

a(x + y) a(y - x) 

p(x) - p(y) = 

a z {x)a z (y) 

or: 

$ K (x)$_ K (x) = p(x) - p(n) 
Furthermore, we recall the following basic addition rules 

C(x + y)- CO) - C(y) = t;^-^ rr > 

2 pix) - p{y) 

and 

p(x) + p(y) + p{x + y) = {C{x + y)- Q{x) - C(y)) 2 = \ ~ P } V ] 

4 V p( x ) - p(y) 



3.2 



3.3 



3.4 



3.5 



3.6 



3.7 



3.9 
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All identities needed in what follows rely on these the above basic addition formulae for the 
Weierstrass family. 

3.2. From addition formulae to PAEs. We will now show that the identites given above 
lead directly to an interpretation in terms of basic elliptic solutions of certain lattice equations. 
Setting 

t = T nym ^ h = a(£) with £ = € + ri6 + in£ + h'y, (3.10) 
and using the abbreviations 

7" — T"n+l,m,h i 7" — 7~ n ,m+l,/i ; 7" 1~n,m,h+l ! 

we can write the three-term relation as: 

a (6 — e)<r(7)r r + a(e — r y)a(5)TT + (7(7 — 5)a(e)rT = , (3-11) 

which is the famous discrete analogue of generalised Toda equation, or Hirota-Miwa equation, 
cf. [2T]. It is probably significant that the fundamental addition formula for the Weierstrass 
c-function obeys naturally one of the fundamental partial difference equations (PAEs) in the 
theory of integrable systems. We will demonstrate now that other PAEs, closely related to cases 
of the ABS list, emerge naturally from the addition formulae as well. 
Thus, setting w = w n ^ m = £(£) — n((5) — m((e) we get from 



(w-w)(w -w) = p{5) - p(e) , (3.12) 

which is the lattice potential KdV equation, cf. |2U[ I3(J|. or the HI equation in the ABS list, by 
the identification p{5) — e := p 2 , p(e) — e := q 2 where e = p{ui) being a branch point of the 
Weierstrass curve: 

r = {(x, y)\y 2 = 4(rr -e)(x- e')(x - e") } 

Setting v(a) = v n ^ m (a) = and identifying parameters p a = $> a (5) , q a = 3><*(e) we also 

have 

p a v(a) - q a v(a) 



v(a) 

p a v(a) + q- a v{a) 



w — w (3.13a) 

= w — w (3.13b) 

v(a) 

from which we get that v(a) obeys a quadrilateral equation which generalises the lattice potential 
MKdV equation, namely 

p a v(a)v(a) + q- a v(a)v(a) = q a v(a)v(a) + p- a v(a)v(a) (3-14) 

Furthermore, setting s(a,(3) = $ Q , +( g(^) , we have 

v(a) v((3) = p a s(a, (3) - p-ps(a, 0) , v(a) v((3) = q a s(a, (3) - q-ps(a, (3) (3.15) 

from which we can deduce: 

j a s - p-ps)(pi3S - p- a s) = (q a s - q-ps)(qps - g_ a s) (3.16) 

which by a point transformation can be shown to be equivalent to (Q3)o, i.e. the Q3 equation 
with parameter A = 0. 
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3.3. Elliptic (non-germinating) seed solution for Q4. We have observed that the addi- 
tion formulae for elliptic functions, and combinations thereof, can be naturally interpreted as 
quadrilateral lattice equations of KdV type. Thus, the elliptic functions constitute elementary 
solutions of these lattice equations. Such solutions can be viewed as seed solutions for Backlund 
chains yielding more complicated solutions which could be considered to be elliptic analogues 
of the soliton solutions. In the remainder of the paper we will construct entire families of such 
solutions, not only for lattice equations of KdV type, but for all equations in the ABS list of 
[2] with the exception of Q4 which is treated elsewhere, Here, however we will restrict 

ourselves by showing that also Q4 admits elementary solutions in terms of elliptic functions, 
and we will do that by presenting the so-called three-leg formula for the equation in the form 
(|2.1d[) . This particular form which essentially is the Jacobi form of Q4, first derived in [18] , 
and for which the 3- leg form was given in [T2], admits a natural parametrization in terms of 
Weierstrass functions as follows. 

Noting the periodicity properties of the Weierstrass functions, i.e., 

o-(t + 2u) = -e 2 ^+^o{t) , C(£ + 2u/)=C(0 + 277 , p(£ + 2w) = p(£) , (3.17) 

in which rj = C( w )i an d w one of the halfperiods, we can introduce the function 

W(x) = ^(x^* , (3.18) 

which obeys the relations 

w(0w(?) = w(t; + O[((0 + ((O + v- ((£ + ? + , w(-0 = -w(0- (3.19) 

Furthermore, we have 

e VUJ 

W(£)W(S + u) = —z n - , W\0 = P(0~e, (3.20) 

and hence as a consequence 

(p(£) " e)(p(£ + u) - e) = (e' - e)(e" - e) = g , (3.21) 

where e' = p(u/), e" = p(w+u/) = —e—e' are the other branch points of the standard Weierstrass 
curve y 2 = 4(cc — e)(x — e')(x — e") . The crux of the matter is the following statement 

Proposition 3.1. For arbitrary variables X, Y , Z we have the following identity 

(X - W(£ + 5)) (Y - W{£ - e)) (Z - W{£ - 5 + e)) 

(X - W(t - 5)) (Y - W(£ + e)) (Z - W(£ + 5 - e)) = 

W{5) (Y W{£) + XZ) - W{e) (X W(£) + YZ) 



(3.22) 



in which 



<r(Z-6)a(Z + e)(T(Z + 5-e) ^ W(S-e) 



ELLIPTIC SOLUTIONS OF ABS LATTICE EQUATIONS 



9 



Proof. Eq. (|3.22p can be established by direct computation through identities that hold for the 
coefficients of each of the monomials XYZ, XY, XZ, YZ, X, Y, Z and 1. The corresponding 
coefficients boil down to the following identities: 

W{£ + e-5) W(Q - W{5) W{e) 

W{i + 5-e) W(0 - W(6) W(e) 

W(£ - e) W(Q - W{8 - e) W{8) 

W(£ + e) W(f) - W(5 - e) W{5) 

W{j + e-5) W(j - e) - W(5 - e) W(e) 

W(£ + 5-e) W(C + e) - W{5 - e) W(e) 

W{j + S) W{j - e) + W{8) W{e) 

W($, - 5) W{i + e) + W[5) W(e) 

g + W{j + 8) W(£ - e)W{j + e - 8) W(Q 

g + W(t-8)W(Z + e)W(Z-£ + 8)W(Z) ' 

and which are all equal to the form given in (|3.23p . All these identities can be proven by using 
the basic addition formula (|3.19j) . the only exception being the last equality which follows from 
the following version of the three-term relation for the cr-function: 

cr(£ + 5) cr(£ - e )a(£ + e- 5) a(£) - a(£ - 5) a(£ +e)a(£-e + 5) cr(£) = <r(2f) a (5 - e) a{8) a(e) , 

and a similar formula with £ replaced by £+w, together with the identification g = e 2r,u / a (ui) and 
the periodicity property (|3.17p of the cr-function. □ 

Identifying in (13.22P the parameters 

P = W{8) , Q = W{e) , p = ((8 + u J )-a8)-r l , q = C(e + u) - ((e) - V , 

which implies that the points p = (P,p), q = (Q, q) lie on the elliptic curve 

f = {(P,p) £C 2 \p 2 = P 2 + 3e + gP~ 2 } , 

cf. [28] . and identifying X = W(£), we see that the expression in brackets on the right-hand 
side of (|3.22p can be written in terms of the following quadrilateral expression 



~ ^ ~ ~s ~ P 2 - Q 2 

Qp „(u,u,u,u) := P(uu + uu) — Q(uu+uu)-\ 

p + q 



uu + uu) — — — ( g + uuuu 



(3.24) 



It is not hard to see that the equation Q^Ju, u, u,u) = is, up to some simple scaling transforma- 
tions, equivalent to the Q4 equation in the form (|2. ld|) . As a direct corollary of Proposition lBTTl we 

have now that u = W(£) = W^ + nS + me) , X = u = W(£) ,Y = u = W® ,Z = u = W(£) , 
is an elliptic solution of the Q4 equation, albeit a trivial one in the sense that as a seed solution 
of a Backlund chain it is non- germinating, in the sense of the discussion of [6]. In the latter pa- 
per we have constructed germinating seed solutions for Q4, in a slightly dfferent but equivalent 
parametrisation, as well as 1-soliton solutions. In a recent paper [11] we have extended these to 
A^-soliton solutions for Q4. 
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4. Cauchy Matrix Scheme 

We will develop now a scheme along the lines of the paper [23] for elliptic soliton solutions, 
based on elliptic Cauchy matrices. In this section we derive the basic relations, and in the next 
section we will use these relations to find a general elliptic TV-soliton solution for Q3. 

4.1. Basic ingredients. At this point let us introduce the Lame function 

:= $ 5 ( K )er^) K , (4.1) 

which coincides with f|3.2|) up to an exponential factor, breaking the symmetry between the 
argument of the function and the su The basic identities for the \& function are the following: 

tt £ («)tf,(A) = e %K ^ +< 5( K )* (5 (A- K )+e %A ^( K -A)^ +( 5(A) , (4.2a) 

* £ («)* 5 («) = e^s +s (K) [£(0 + CO*) +C(«)-C(£ + * + *)] , ( 4 - 2 b) 

* C («)^(A) = ^(« + A) [C(0 + C(«) + C(A)-CK + « + A)] , (4.2c) 
in which we have introduced 

« = ««) = C« + Q - C(«) - C(fl = 1 • (4-3) 

Furthermore, we have the symmetry: ^?s(— n) = —^>s(k) . 

The starting point for our construction is the "bare" non-autonomous Cauchy matrix 

M° = (M^) . . =1) jJV , vKiK, • k,). (4.4) 

depending on a variable £ which is linear in the independent variables n, m, namely £ = £o + 
nJ + me , with <5, e being the corresponding lattice parameters. We will assume that the set 
rapidity parameters {/tj, i = 1, . . . , iV} is such that + Kj 7^ (modulo the period lattice of the 
Weierstrass functions) . 

Furthermore, we redefine henceforth the lattice parameters p K , q K of section 3, to include an 
exponential factor, and thus we define 

p K = <& s (k) , q K = * e («) (4.5) 

Setting k = i/Cj, k = ±/tj, we can derive from the basic addition formula (|4.2aj) the following 
dynamical properties of the elliptic Cauchy matrix: 

p Kl Mlj = m 5 {ni)^^i + Kj) 

Similarly, we have 

P- K Ml 3 = M° jPK] e-^+^ - e -^^^( Ki )* e+5 ( Kj ) . 

1 Although most of the results of this paper can be obtained in terms of the ^-function alone, working with 
the function \1/ has certain advantages as the latter is the more natural function in connection with associated 
continuum equations. We note, however, that the inclusion of the exponential factor amounts to a specific gauge 
transformation on the quantities defined later on in the constructions, and hence could be removed without 
affecting the main results. 
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We introduce now the plane- wave factors (i.e., discrete exponential functions) 

' a -rirX' f4^4)"*<= 2 ««"Po,oW, (4.6) 

and for the specific values Ki of k the quantities pi := p(ftj) obeying the shift relations 

El — c 2y S Ki P-Ki ^ Pi_ _ e 2T) e m 9-Kj ^ / 4 

Pi p Ki ' Pi q Ki 

where the superscripts ~, ~ as before denote the elementary lattice shifts with regard to the 
independent variables n and m respectively. Use has been made of (|4.3|) . noting that £ = £ + 8 , 
and a similar relation for n e using £ = £ + e . 
Next we introduce the A^-component vectors 

»• = (Pi*e(«i)) i= i,...,jv . s = ( c j*e(«j))j=i,...,jv ( 4 - 8 ) 

where Cj are iV constants w.r.t. the discrete variables n, m, and pi = p n ,m( K i) , in terms of 
which we define now the "dressed" Cauchy matrix: 

M = !A7,.,)„, , v , M itj = piM)' jCj . (4.9) 

As a consequence of the relations given earlier, and employing the definitions of the plane- wave 
factors, we can now describe the discrete dynamics as follows: 



Lemma 4.1. The dressed Cauchy matrix M, defined in (|4.9p . obeys the following linear rela- 
tions under elementary shifts of the independent variables n 

e T,t > K p^ K M - Mp^ K e VsK = rs T (4.10a) 
Mp K e^ sK - e~^ K p K M = rs T , (4.10b) 

and under shifts of the variable m the similar relations: 

e^ K q^ K M - Mq_ K e^ K = rs T (4.11a) 
Mq K e~ VsK - e~ VeK q K M = rs T . (4.11b) 

Here we have adopted the matrix notation, with K = diag(«i, . . . , kn) denoting the diagonal 
matrix of the Ki parameters. Similarly, the symbols p±k, Q±k denote the diagonal matrices 
with entries p± Ki , q± Ki . 

In what follows we will employ the relations (|4.10p and (|4.1ip to obtain nonlinear discrete 
equations for specific objects that we will subsequently define in terms of the Cauchy matrix 
M. 



4.2. The r-function and related basic objects. Introduce now the r-function : 

/ = fn,m = det (1 + M) , (4.12) 
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where 1 is the N x N unit matrix. From the relations (j4. lOj) and (|4.1ip we can derive the 
discrete dynamics of the r-function by performing the following computation 

/ = det (l + M\= det {l + [e™ K p_ K M - fs T ] e~ %K (p-k)' 1 } 

= det {e^ K p^ K [1 + M - e-™ K (p- K )- l rs T ] e~ raK {p- K )- x } 
= det { (1 + M) [1 - (1 + M)- 1 e-™ K {p_ K )- 1 rs T ] } 
= / det{l- {l + M)~ 1 e~^ K {p- K )- 1 rs T } , 

and similar computations involving the other lattice shift. Noting that from (14. 7\i it follows that 

e-^ip^r = e ^^|i^| r , (4.13) 

we can use (|4.2c|) and the Weinstein-Aronszajn formula 

det (l + a b T ) = 1 + b T ■ a 

to write 

f = 1 - a T(i + M)" 1 [C(0 + C(K) + C(S) -((Z + K + S)}- 1 r = V 5 . (4.14) 

Here and in what follows the notation C(-K") denotes the diagonal matrix with entries C( K i); 
(i = 1, . . . , iV), whilst all scalar terms are supposed to be considered as coefficients of the unit 
matrix. In (|4.14j) we have introduce the quantity: 

V a = 1 - s T (1 + M)" 1 ( Xq , K Y l r = l-s T ( Xa , K + M)- X r , (4.15a) 

for arbitrary parameter a. In addition to V a we will also need, for later purpose, the objects 

U a ,{) = * T (Xf3, K ) _1 (1 + M)- 1 ( Xq , K )~ X r = U^ a , (4.15b) 

where both a, are arbitrary (complex) parameters. In (14.15b|) and (j4.15aj) we have abbreviated 

Xa,p = XaAO ■= C(o) + C(P) + C(£) - C(£ + a + P) , (4.16) 
and where Xa,K denotes the diagonal matrix with entries Xam- 

4.3. Basic linear relations. In order to derive relations for the objects V a and U a ^ we need 
to introduce the iV-component column- resp. row vectors: 

u a = (l + M)- 1 ( X a,^) _1 r , (4.17a) 
% = s T (xp, K )- 1 (l + M)- 1 , (4.17b) 
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Performing the following calculation: 

u a = (1 + M)- 1 ^,*)- 1 ? = (1 + M)-\xc l ,K)- l e 2 ^ K \f K ^ K r 



e 



(1 + M)p K e" nsK - rs T 



OK) + C(g + a) - C(g) ~ C(K + e + a) 

C(*0 + C(£) + C(«)-C(^ + £ + a) 



f 1 C(<*) + C(Q-C(g + <*) + Cffl ) 

" ^ C (K) + c(e) + c(«)-c(^ + e + a)J 

multiplying both sides by (1 + TVf) -1 and introducing the vector 

u° := (l + M) -1 r , (4.18) 

we get the relation 

e~ risK p K u a = (s T u a )u° + Xa,su a -u° e~ r ' sK p K u a = -V a u° + Xa,8U a . (4.19a) 
In a similar way one can derive the relation: 

e^ K p^ K u a = V a u° - Xa,-su a . (4.19b) 

A similar set of relations can be derived for the adjoint vectors ()4.17b[) which involves the adjoint 
vector to (|4.18p . namely 

V := s T (1 + M)' 1 , (4.20) 

and obviously these relations all have their counterparts involving the other lattice shift related 
to shifts in the discrete independent variable m instead of n. 

Summarising the results of these derivations, we have the following statement: 

Lemma 4.2. The N -component vectors given in (|4.17p . together with the ones defined in (|4.18p 
and (|4.20p obey the following set of linear difference equations 

e~ nsK p K u a = -V a u + Xa,su a , (4.21a) 

e^^-p-KUa = V a u - Xa,-5U a , (4.21b) 

^^xe^ %K = Vp l u - x/3-S (4.21c) 

% p- K e mK = -Vp *w + X p,S 'up , (4.21d) 



and a similar set of relations involving the shifts in the variable m obtained by replacing p±K 
by q±K, S by e and ~ by ~. 
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4.4. Basic nonlinear relations. From the basic relations (|4.2ip we observe that 

_^ T ^ + p{K)^ 5 {K) g _ gT C(jQ + C(*) + C(g + g) ~ ffl + g + g) ~ 
^K)^(X) * C (K) + C(/ 3) + C(e) _ C(e + K + /3) " c 



-T 



1 + — — — -~ -~ u a = (1-V a ) - X0,-sU a ,p 



= Xa,sU a ,l3 -V a (l- Vp) , 
from which we get the following relation: 

V a Vp = l- Xp,-sU a ,p ~ Xa,sU a ,p . (4.22a) 

Similarly we have 

V a Vp = l- Xp,-eU a ,p - Xa,eU a ,p . (4.22b) 

These relations can also be rewritten as follows 

V a ) (^(/3) Vp) = Pa e-^ a [*t(a + /?)(1 - X<x,pU a ,p)]- P -p^ [%(a + - Xa,/?^) 

(4.23) 

provided a + (3 ^ (modulo the period lattice of the Weierstrass functions) . 

We now intend to remove the ^-dependence in the coefficients of (|4.23p by performing yet 
another change of variables, namely by defining 

W a := ^(a) V a , S a> p :=^^(a+(3)(l-Xa,pU a ,p) = Sp ia , a+(3 / (mod. root lattice, 

(4.24) 

We now have the following result: 

Lemma 4.3. The following difference relations hold between the quantities defined in (|4.24p 

W a Wp = p a e-^ a S a ,p - p.pe^S a ,p , (4.25a) 
W* Wp = q a e-^ a S a ,p - q-pe^S a ,p , (4.25b) 

where and a + (3 ^ (mod. period lattice). Furthermore, in the latter case the relations (|4.25p 
are replaced by the following ones: 

W a W. a = p a e~i* a {[C(0 - C(9 " C(«) + C(« + ^)] + - p(a)}U a ,. a - [p(0 - p(a)]U a ,- a } , 

(4.26a) 

W a W. a = q a e-^ a [[((0 - C{0 ~ ((a) + C(« + e)] + [p(0 - p(a)]U a ,. a - [p(0 - p(a)]U a ,. a } . 

(4.26b) 

As a consequence of the relation (|4.26p . taking the limit a — > —5 and a — > —5 respectively, 
we obtain the inversion relations 

~ pK(5-c(0]« P [c(o-c(0]e 
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4.5. NQC equation. We now present the nonlinear PAE for the variables S ai p which follow 
directly from the relations (j4.26f) . in fact by simple elimination of the W variables. 

Proposition 4.4. The variable S a ^, defined in (j4.24j) . for fixed parameters a, (3 obeys the 
following nonlinear partial difference equation 

(p a e~^ a S a ,p - p-pe^S a ,p) (ppe-^S a ,p - P - a e^ a S a ,p) 

= (q a e~^ a S a ^ - q.pe^S^p) (qpe'^S^p - q. a e^ a S a ^) . (4.28) 
Furthermore, the following relation holds 

= (qfye-TP Sw - q^%^) {q^e^S^ - q- a e* a X,a>) (4.29) 

between these solutions of the partial difference equation (|4.28j) with different fixed parameters 
a, /?, a' , /?'. 

Eq. (|4.28|) first appeared in a slightly different form in [30| . It should be noted that the latter 
parameters a,(5,a',f3' are to be distinguished from the lattice parameters 8 and e which are 
associated with the lattice shifts. 

Eq. (|4.28p can be cast in a more universal form by setting u a $ = p 1 ^ 2 (a)p 1 / 2 ((3)S a ^ , which 
leads to an equation of the form: 

P(u a ,i3Ua,i3 + u aj pu aj p) - Q(u a> pu a ,p + u a ^u a ^) = (p 2 - q 2 )(u a ^u at p + u a> pu a ,p) , (4.30) 

i.e. the Q3 equation for A = 0. Here we have used the identities p a p~ a — q a Q-a = PpP-p — 
qpq~p = p 2 — q 2 = p{5) — p(e) and we have introduced the parameters P, Q obeying 

P 2 = PaPpP-aP-p = {p 2 ~ a 2 )(p 2 - b 2 ) , (4.31a) 
Q 2 = q a qpq- a q^ = (q 2 - a 2 )(q 2 - b 2 ) , (4.31b) 

i.e. p = (p, P), q = (q, Q) are points on a (Jacobi) elliptic curve, with moduli ±a, ±6 which are 
given by a 2 = p{a) — e , b 2 = p{(3) — e . From (|4.29j) it follows also that any four solutions u a> p, 
V,/?', ««,a', u p,f3' are connected via the relation 

Pa,a'Ua,pUa' ,/3' ~ P<y! ,pU a ,pU a ' ,/?' ~ Pa,/3'U a ,pU a > t p> + Ppfi>U a ^U a i^i = 

= Qa,pUa,a'Up,p> - Qa' ,pU a ,a'Up,P' ~ Qa,P'U a ,a'Up,p' + Q a > ,p'U a ^ a iUpfii , (4.32) 

in which the parameters P ai ,a 2 i etc are lattice parameters associated with different elliptic 
curves, i.e. they obey relations of the type 

PL a2 = {p 2 -a\){p 2 -a 2 ) , Ql ua2 = (q 2 -a 2 )(q 2 -a 2 ) 

with branch points dbai, ±02 which are associated with uniformising parameters ax, 02 (taken 
in the set {a, a', (3, /?'} ). The various lattice parameters are not independent, but are related 
through the quadric relations 

P 2 , a > ~ Pl>,p ~ Pl,p> + Pp,p> = (« 2 - b 2 )(a' 2 - b' 2 ) , 
and a similar relation for the parameters Q. 
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4.6. Miura relations. In addition to the relations between the objects W a and S a ^ we need 
another set of relations involving the object 

w := C(0 - n((6) - mCO) - s T u , (4.33) 

which, as we shall see, solves the HI equation. 

Lemma 4.5. Between the objects W a and w the following relations hold: 

p a e-^ a W a - q a e~^ a W a P - a e^ a W a - q- a e^ a W a 



w — w 



w — w 



p a e~^ a W a + q_ a e^ a W a p_ a e^ a W a + q a e~^ a W a 



W a W a 



and, furthermore, w obeys the equation HI in the following form: 

(w — w)(w — w) = p 2 — q 2 , p 2 = p(S) — e, q 2 = p(e) — e 
Proof. From the relation 



(4.34a) 
(4.34b) 
(4.35) 



e VsK pKU a = Xa,SU a - V a u 



e ''""PKUa = s T J) T J^ 5 (K) e - VsK u c 



Xa,5S T U a - V a S T U 



we get 



C(K) + C(£) + C(S) - ((K + u a = (l- V a ) Xa ,s - w V a , 



where wq = s t uq . Similarly, from p-KU a = V uq — Xa,-8 u a we can derive the relation: 

C(K) + C(0 - t(8) - C(K + £)] u a = V a w - (1 - V a )xa,~s ■ 
Combining these with the analogous relations involving the other lattice shift, i.e. 

C(K) + C(0 + ((e) - ((K + 0] u a = (1 - V a ) Xa ,e - w V a , 
C(K) + c(0 - COO - C(K + 0] u a = V a w - (1 - V a )x a ,s , 



to eliminate the terms involving K on the left-hand sides, we obtain the following set of relations: 



w — w 



w — w 



Xa,8Va ~ Xa,eV a Xa-eVa ~ Xa,sV a 



Xa,sV a ~ Xa-e V a _ Xa, £ V a - Xa,-S V c 



(4.36a) 



(4.36b) 



V a V a 

The relations (|4.34|) are obtained from (|4.36|) by substitution for V a using (|4.24j) . and by using 
the expressions for Xa,±6 an d Xa±e i n terms of the \F-function in the form 

^p- a ^^a) 



A«,0 c , T , / \ J Act, — o ° 



* ? (a) 



*f(a) 
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and similar expressions with 5 replaced by e. The latter relations follow directly from the 
addition formulae (|4.2j) . Finally, setting a = 5 p~ a = p~$ = , the second forms of the 
expressions in the lemma yield: 

w — w = —Q-ge ' 77— and w — w = q$e le -^=— 

which by multiplication, using q$q-s = p(e) — p(5) gives us the HI equation for w. □ 

5. Elliptic A^-soliton solution of Q3 

5.1. iV-soliton formula for Q3. We are now ready to formulate the main result of this paper, 
which comprises a generalization of the rational iV-soliton solution for Q3 as obtained in [6j [23] 
to the elliptic case. We use the following expression to abbreviate the quadrilateral form of the 
equation: 

Q^ q (u, u, u, u) := P{uu + uu) — Q(uu + uu) — {p 2 — q 2 ) f(uu + uu) + ~pQ j (5-1) 

in which the lattice parameters p, q are points (p,P), (q,Q) on the elliptic curve: 

r Q , b = {(x, X) I X 2 = (x 2 - a 2 )(x 2 - b 2 ) } , (5.2) 
with moduli ±a, ±6, which we refer to as the parameter curve. 
Theorem 5.1. The following formula 

«w = V /2 («)p 1/2 (/?)s Qi/3 + V /2 («)p 1/2 (-/3)s q ,- /3 

+Cp 1 / 2 (-a)p 1 / 2 ((3)S- a , f 3 + Dp 1 / 2 (-a)p 1 / 2 (-P)S- a ^ , (5.3) 

for each fixed positive integer N , and parameters and variables given in the previous section, 
provide solutions of the quadrilateral equation 

Q^(u( N \d N \rt N \^ N) )=0 (5.4) 
with arbitrary constant coefficients A,B,C ,D, where A is given by 

A = p'(a) p'ip) det ( £ ■ (5-5) 

Furthermore, the associated function to : 

rj(N) = Ap 1 / 2 (a)p 1 / 2 ((3)W a W^ + Bp 1 / 2 (a)p 1 / 2 (-P)W a W_ f3 

+Cp 1 l 2 {-a)p 1 ' 2 {0)W- ol Wp + Dp 1 l 2 {-a) P x ' 2 {-0)W- ol W-p , (5.6) 
factorises the corresponding biquadratic to the quadrilateral, namely it obeys 

X 9 (uWtf*r)) = U Wum , M q (#),2W) = C/( JV )^( JV ) (5.7) 

where 

Oi p (u, u) := P(u 2 + u 2 ) - (2p 2 -a 2 - b 2 )uu + — (5.8a) 

?Cq(u,u) = Q{u 2 + u 2 ) - (2q 2 - a 2 -b 2 )uu + ^- . (5.8b) 
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Proof: The proof is obtained through direct computation, using the relations between the various 
objects established earlier. It follows roughly the same lines of computation as in the rational 
case [24] , but only differs in a few subtleties regarding the elliptic representation. It breaks down 
into three steps. 

step # 1: Using the basic relations established between the objects W a , Wp and S a ^ one can 
obtain relations between the objects ti™, and the HI object w, as given in the following 

lemma. These constitute what in effect is a Miura transformation betweem Q3 and HI. 



Lemma 5.2. The following relations hold between the elliptic N-soliton solution u 



Q3, together with its associated object U = U^ N \ and the elliptic N-soliton solution w of HI: 

Pu - Qu - (p 2 - q 2 )u Pu — Qu- (p 2 - q 2 )u 



of 



w — w 



w — w 



U 



U 



Pu — Qu — (p 2 — q 2 )u Pu — Qu — (p 2 — q 2 )u 



U 



U 



(5.9a) 



(5.9b) 



Proof. These relations, which can be obtained by direct computation term by term in the ex- 
pressions (15.3j) and (15. 6ft (after multiplying out the denominators), namely as follows: 



Pu — Qu — (p 2 — q 2 )u 



Af /2 (a)f /2 ((3) 



p[ e 2((8)aP^\ ( e 2C(5)/3_PM e hC(5+<5)K(0] £ 
P-aJ V P-p) 



Q e 



1/2 



3 2C( £ )/3_^_ 

q-p 



1/2 



Ap {a)p (J3) 
Af /2 (a)f /2 (f3) 
Af /2 (a)f /2 (f3) 



P, 



tPpe °a,p 



e i-C(f+e)+C(€)](a+/»)5 a . 



(p 2 - q 2 )S aj p 



ppe 



-vsP 



p a e 



-■q s ag 



a,/3 -P-pe 



VsPcj 



pp 



e -vsP W W 



a Wp- qpe 



-rfeP 



W a Wp 



a,P 



+ 



Qpe 



-VeP 



(p 2 - q 2 )S a ^ 

+ •■■ 

P -Veaq 



+ 



q-p( 



Ap (a)p (P)(w-w)W a W (, + ■■■ 



by using in the last steps the basic relation (14.251) in combination with (|4.34j) . and where the 
• • • stand for similar expressions for the remaining terms with coefficients B, C, D instead of A, 
and with (a, (3) replaced by (a,—f3), (—a, (3) and (—a,—f3) respectively. Thus, identifying the 

right hand side, up to a common factor (w — w) with the variable U, we obtain the first relation 
in (|5.9aj) . The other relations in (|5.9j) follow by similar computations. □ 



step # 2: We next establish the main identities between u = u^m, and U = lli N m. 
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Lemma 5.3. the following factorisation formulae for the biquadratics associated with the Q3 
quadrilateral hold: 

UU _ P(u 2 + ~2 } + {2p 2 _ a 2_ b 2 )u ~ = P'(<x)p'(P) det{A) ^ (5 1Qa) 

UU -Q(u 2 + u 2 ) + {2q 2 -a 2 - b 2 )uu = ti^M. det (A) , (5.10b) 
in which the 2x2 matrix A is given by 

A= ( n n ) ■ ( 5 - n ) 



v C D , 

Proof. This can asserted by computing the (Lax type) matrices 

_ ( Pu-{p 2 -b 2 )u , e^p.pU \ 

V e-nPppU , -Pu + (p 2 -b 2 )u J ' [bAZa) 

M - V e-^PqpU , -Qu + {q 2 -b 2 )u ) ■ {b ' Ub) 

in two different ways: one directly, and the other by expressing the entries in terms of our basic 
objects. In fact, introducing the 2-component vectors: 

rl = ( p V 2 (a)W a ,pV 2 (-a)W_ a ) , r p = ( //^ ^ ) , (5-13) 
we can express L as: 

L= ( -^WP^rlArp , ei'Pp-prlArp \ , . 

V e-^pprlArp , -^pppZpf^Arp ) ' 1 ' ' 

and a similar expression for M. In fact, 
Pu - (p 2 - b 2 )u = 
= V /2 («)p 1/2 (/5) [p-aP-^^S^ - PffP-pScf,] + ... 
= ApW^pWipW'tp-p (p-ae^S^ - ppe^ScJ) + . . . 
= _ V /2 («)p 1/2 (/?)e r « V/^W, + . . . 

= -v^^ (Ap l ' 2 {a)p l l 2 (f3)W a Wp + •••)= -y/p^pIprlArp , 

and similarly, 

Pu - (p 2 -b 2 )u = yjppp-p r^Arp . 
Evaluating the determinant of L using the general identity 

det f x jVj\ = det {ivl ■ x j)>:i I ) 
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holding for any 2r r-component vectors x, y, we now get: 
det(.L) = P/3P-/3 det ^Arp r T a + Arp r 

= {p 2 - b 2 ) det (A) det (r p r T a + r p r T a 



(p 2 -tf) det{A) det i(rp,r,3) 



T 
a 
-T 



= -(p 2 - 6 2 )det(7l) det (r a ,r a ) det (rp,rp) . 

It remains to compute the determinant of the matrix (r a ,r a ) whose columns are the 2- 
component vectors r a and r a . This is done by using (|4.26p as follows: 

det (r(a),r(a)) = p 1/2 {a)p 1/2 {-a)W a W- a - p 1/2 {a)p 1/2 {-a)W a W. a 



Pa 



e~ Vsa W a W- a 



P-a 



P~a V Va 

= ^p-Z^[2C(a) + C(5-a)-C(5 + a)} 
Thus, putting everything together we obtain the result: 



P'(a) 



det(-L) = -{p 2 -b 2 )det(A) 



Vp 2 



\Jp 2 — b 2 



On the other hand a direct computation of the determinant gives: 

det(L) = -[Pu- {p 2 - b 2 )u] [Pu - {p 2 - b 2 )u] - {p 2 



b 2 )UU 



(P 2 



P(u 2 + u 2 ) 



(2p 2 -a 2 - b 2 )uu - UU 



which yields the first of (|5.7p . the second part of which follows by a similar computation of 
det(iVf), which is obviously obtained from the det(L) by replacing 5 by e and ~ by ^. □ 

step # 3: The final step is to arrive at Q3 by combining the results of step # 1 and # 2 for 
the solution from the Miura relations (|5.9p . together with the relations (|5.7p . Thus, since 
w obeys the HI equation, we have: 



(w — w)(w — w) 

Pu — Qu — (p 2 — q 2 )u Pu — Qu — (p 2 — q 2 )u 



1 

UU 

2 



(P 



Q(u 2 +u)- (2q 2 



b 2 )uu 



Q 



P l (uu + uu) + (Q z - P z )(uu) + (p z - q 2 yuu - PQ(uu + uu) + 
~{P 2 - q 2 )P(uu + uu) + (p 2 - q 2 )Q{u 2 + u ) 



which after some cancelations leads to Q3 in the form (|5.4p and, hence, concludes the proof of 
the Theorem. 
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5.2. Hirota form of the elliptic iV-soliton solution. The r-function can be explicitly com- 
puted by using the expansion 



N 

/ = det(l + M) = l + £|M iii | + £ 



i=l i<j 



H hdet(M) . 



and using the Frobenius formula for the relevant elliptic Cauchy determinants, cf. [16]. Thus, 
from 

-*<*w*> - (n^j) *^«-™^n(^)° ■ (-> 

Introducing the notations 

e A itj ._ ( ~ K i) \ 2 ; e e x = Pi* C -2C(£U,: ^ 

the Hirota formula for the r-function thus takes the form: 

/ = i+E^I^W E ^+ 2 g+ 2 ^) e ^ + A, + (516) 

i=l CJ ^'' «<j=l ^ ' 



A? 



E cr(C + + 2Kj + 2K k ) j. +e . +ek+A , +A . k+ A jtk , 
i<J<fc=l vsy 



(5.17) 



The expression ()5. 17|) for r-function / enters solution of Q3 through the main quantities S a r 
and W a in the following way. First, we note that V$, and hence W$ can be expressed in terms 
of / via (|4.14p . Second, by setting a = e, j3 = 5 in (|4.22ap . using also (|4.14j) . we obtain 

S 5 , E = Vz(8 + e)t . (5.18) 

Thus for special parameters a, (5 associated with the lattice parameters, we have an explicit 
expression for these quantities in terms of the r-function. To get a similar expression for arbitrary 
a, (3 we just need to extend the two-dimensional lattice in terms of the variables n and m 
(associated with lattice parameters 5 and e respectively) to a four-dimensional one, containing 
additional lattice directions associated with the parameters a and (3. The discrete variables 
associated with these parameters we will denote by h and I respectively, and hence the extension 
amounts to including in the plane-wave factors p(n) a dependence on these variables as follows: 



(5.19) 

whilst extending at the same time the variable £ covariantly as £ = £o + n5 + me + ha + 1(3 . In 
(|5.19p we have set a K = ^ a {n) and b K = ^(/c). In this extended lattice the elementary lattice 
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shifts T a , Tn associated with these new lattice directions can be used to express S a ,/3 and W a in 
terms of the r- function, in a way similar to (15.18P and (|4.14j) . leading to the expressions 

S a ,p = ^(a + (3)^1- , W a = ^(a)^ , (5.20) 

and where the T a , Tp denote elementary lattice shifts in the directions associated with the 
parameters a, (3. 

Combining now the expression (|5.20[) with the explicit form of the the solution (|5.3p and of 
its associated variable (|5.6p we obtain the following result. 



Theorem 5.4. The elliptic N-soliton solution (|5.3p of Q3 can be written as w N ' = (7 A f) / f 
where the shift operator 7 A is given by 



J A f .- Lr^r p (r a ) | c ^_ a + p)T -i Tfj , D ^- a _ f3)T~ l T~ l ^> 



whereas the associated solution (|5.6p can be written as = (T>j±f • /) / f 2 where the bilinear 
operator D A is given by 



V A f-g := m 



- LI S'/3l'aJ I /"l,Tf I „AiTr .//S\T-lo,T. n,T, / „,\iTr . / flAT"— 1 o T^-l I f w 



'(5'22) 

luii/i m denoting the pointwise product in the space of functions of the lattice sites, making use 
of the notation (I5.13p . (I5.1ip as well as 



(5.23) 



5.3. Linear scheme for elliptic soliton solutions. We will now show how the corresponding 
solutions constructed in section 4 arise naturally from a scheme of linear equations, which in 
turn leads to a novel Lax representation for Q3. These also explain the emergence of the Lax 
type matrices that were used in the proof of Theorem 15.11 

The objects (|4.17p form the basis of the linear structures underlying the lattice equations that 
we have studied in the previous sections. In fact, they form the main ingredients, on the basis 
of which we can derive Lax pairs for those lattice equation^. 

Associated with the Cauchy matrix M we introduce the following parameter-dependent 2N- 
component eigenfunctions; 

W(a) := (,T(p 1/2 («),P 1/2 (-«))(^ ( V K) VA-l + K))^**)- 1 

(5.24b) 



2 An alternative approach being the one using the multidimensional consistency property, cf. [231 1131 132] to 
derive Lax pairs from the equations themselves, but the resulting Lax pairs often require awkward prefactors 
containing square roots of the relevant biquadratics. 
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where the vectors r° and t s° are given by 

(r°) T = ( Pl ,... , PN ) , ( s °f = ( Cl) ... , CN ) . 

Note that these vectors (|5.24p effectively live in the direct product space of the 2-component 
vectors labelled by the parameters ±a, ±/3 and ^-component vectors associated with the soliton 
parameters Ki, (i = 1, . . . , N). The dynamical properties of these vectors are described by the 
following set of relations: 

Proposition 5.5. The following set of relations describe the behaviour of the vectors given in 
(|5.24|) under the elementary shift in the coordinate n of the multidimensional lattice 

e^ K p_ K u K {(3) = ^p£pZ£u K (P) +r (3 u° K , (5.25a) 
e-^ K PK u K (f3) = ^p^u K {f3)-r p u° K , (5.25b) 

W(a)p_jce' teJC = v^pl^ W(a) - *«K > (5.25c) 
l u K (a)p K e~^ K = ^p—'uKia) + S&r^ . (5.25d) 

where the 2-component vectors and rp are given in (I5.13p . and where the N- component 
vectors u® K and *tt^ are given by 

u° K = (l + M)- 1 r , *u° K = s T (l + M)" 1 , (5.26) 

recalling the vectors r and s as defined in (|4.8|) . TTie shifts in the discrete variable m (and other 
lattice directions) are described by similar relations like (|5.25[> . oy replacing 5 is replaced by e 
(i.e. p±K replaced by q±K), and~ by^. 

Proof. The proof of these relations follow similar derivations as those in section 3, i.e. making 
use of the basic relations (|4.1ip (and (|4.10p ). for the elliptic Cauchy matrix, as well as the 
definitions of the various quantities such the plane- wave factors p, (|4.6|) and the objects (|4.15a|) . 
For instance, to prove (|5.25ap we proceed as follows: 

(1 + M)e™ K p_ K u K {P) = [e™ K p_ K {l + M) -rs T ] u K (/3) 

in which the entries in the matrix term on the right hand side are computed as follows 
enXp-KVt&P + K)p 1 / 2 (±(3)r° = e r ' sK p- K ^ ? (±/3 + K)p 1 ^(±P) e - !b " K —f p 



P-K 



/ \ 1/2 

e -^ K PK ^J±p + K)p 1 / 2 {±(3)f = e^ K y 5 (K)y £ (±(3 + K)p l ' 2 {±3)e™P (EM) jfi 



e -V S K 



e"° K *s + 6(K)*t(±l3) + e ^ ±/3+ ^ 5 (T/?)^ +5 (±/? + K) p l ' 2 {±6)e^ V -±t 

where use has been made of the addition formula (I4.2ah and the shift relations for the plane 
wave factors p(±/3) and in the components of the vector r°. Inserting these expressions in 
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the former relation, and multiplying both sides by (1 + M) 1 , whilst taking note of the fact 
that ty£+s(K)r° = r according to the definitions, we obtain (|5.25ap . by identifying 



s T u K {(3) 



K 



/(1 + M)^WJ 

p l/2 { _ md _ m _ V _^ 







and rewriting these 2- vectors in terms of W±p, using the definitions (|4.24p and f|5. 13j) . Similar 
computations yield the other relations in (|5.24p . where in the latter two we need to indentify 
the combination 



W(a) r = (pV 2 {ct)*t:{a){\ - V a ) , pV2(_ a )^(_ a )(i _ y_ Q )) . 



□ 



5.4. Lax representation. Starting from the basic relations (15.250 for the 2N-component vec- 
tors (|5.24|) we will now define first iV-component vectors by projecting in the 2-component vector 
space by making use of the arbitrary coefficient matrix A of (|5.1ip . Thus, multiplying (I5.25ap 
and (|5.25bp from the left by the 2-component row vector r^A we obtain respectively 



e^ K p. K r T a Au K ((3) = ^pppIprlAu K (fi) + U^vP K 
e~ VsK p K rlAu K (P) = y/ppp-p rlAu K (P) - ^ppP-p u 



(5.27a) 

Vp^^W] <(5.27b) 

where we have used the identities relating the expressions in the entries of the matrix (|5.12ap 
and (|5.14p . Similarly, by multiplying (|5.25ap and (|5.25bp from the left by the 2-component row 
vector r a A we obtain respectively 

e^ K p. K r T a Au K {l3) = ^p—pr T a Au K {(3) + [y/p^u^ - ^p—pu^] w^(5.27c) 

e-^ K p K rlAu K (P) = ^pp-p^7lAu K ((3) - U^u° K . (5.27d) 

Introducing now the new object 

u K (a,j3) =rlAu K {P) . (5.28) 

which is a iV-component vector, and eliminating the quantity r^AujciP) from the pair of 
relations (|5.27a|) and (I5.27bl) we obtain 



(b 2 - k 2 )u K (a,(3) = p K e-^ K U^u° K + [Pu^ - (p 2 - b 2 )u 
whilst from (|5.27cj) and (|5.27dp we get 

(b 2 - k 2 )u K (a,(3) = -p_ K e^ K U^u Q K + fp«W - (p 2 - b 2 )u 



(TV) 



-AN) 



71° 

u K 



u° K 



(5.29a) 



(5.29b) 



These relations can be used to constitute one part of a Lax pair, whilst the other part follows by 
replacing 5 by e (i.e. p±k replaced by q±K), and ~ by ~. Thus, in terms of the 2N-component 
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vector 



u 



K 



uk{ol,0) 

the following new Lax representation for Q3 is obtained 

.PuW + (p2 - &>W 



q K e-^ K 4> K (a,f3) 



b2 -P 2 _ A _ tj(N)tj(N) 
yTZj? P u u 



b 2 -k 2 

Pu (N) _ (p 2 _ 62) ~(7V) 



(TV) 



(5.30) 

t> K (a,/3) 
(5.31a) 
t> K (a,P) 
(5.31b) 

In a similar way we can use the relationns (|5.25cj) and (|5.25d|) to derive an "adjoint" Lax pair 
for the row vectors defined by 

l 4> K {a,l3) = (*u° K , W(«,/?)) , W(«,/3) = W^-Ar^ , 

for the entries of which one can derive the following coupled set of equations: 

W(a,/9)(o 2 - k 2 ) = *vP K p- K &> K uW + *< \Pu^ - (p2 _ a Z) u {N) 



QuW + (q 2 - b 2 )u 
^A_ u (N)fj(N) 



Qu 



b 2 -k 2 

_ ( ? 2 _ 6 2 )2 (7V) 



(5.32) 



(5.33a) 



u if 



W(a,/3)(a 2 -fc 2 
In terms of (|5.30p one can derive a similar Lax pair, namely 

1 / + (r? - n 2 )^ 1 



Pu (N) _ (p2 _ a 2 } ~(iV)l _ t u KpKe -r, s Kfj(N) (g 33^ 



4> K (a,(5)p- K e 



% K (a,f3)q. K e^ K 



>K 



-PuW + (p 2 - 
a 2 -k 2 



-QuW + (q 2 



a?)u 



(N) 



k l 



, p>)_ (j ,2_ a 2) 5 (W) J 

(5.34a) 



(5.34b) 

In (|5.3ip the structure derived from the elliptic iV-soliton solutions, following the computations 
in this section, leads to 2Nx2N Lax matrices with an NxN block structure, each block of which 
is labelled by the soliton parameters Kj. However, having derived the Lax pair, one can now 
consider the variables k (which were diagonal matrices) as general spectral parameters, and 
assess the compatibility of the system (|5.31[) on its own merit. This leads to the following 
statement: 



Theorem 5.6. The compatibility conditions for the two linear equations (|5.31ap and f|5.31b|) are 
satisfied if and only if obeys Q3 in the form (|5.4p . making use of the biquadratic identities 
(|5.7p together with (|5.8p . The same holds true for the two linear equations (|5.34ap and (]5.34bp . 

Proof. The proof is by direct computation, where the (1,1) and (2,2) entries yield the equalities 
at the left-hand sides of (|5.9p . whilst the (1,2) and (2,1) entries yield the equation inserting the 
explicit expressions for the biquadratics "Kp and "Kq of (|5.8p . □ 
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6. Backlund scheme for the ELLIPTIC SOLITON SOLUTIONS 

In [23] we showed by explicit computation that the recursive structure of the iV-soliton solu- 
tions for Q3 coincides with the action of the Backlund transformation. Although this is to be 
expected, the mechanism of the proof reveals the role of the relations between the various quan- 
tities obeying different equations in the ABS list, in particular the Miura type transformations 
between them. In Theorem 16. II we make the statement as in [24], but now for the elliptic soliton 
solutions, and uncover the various aspects of the proof which have to be adapted for the elliptic 
case. 



6.1. Backlund transformation from the N- to the (jV + l)-soliton solution. The fol- 
lowing theorem relates the elliptic (N + l)-soliton solution of the form given in (15. 3ft to the 
corresponding elliptic A^-soliton solution with one less soliton parameter, thus establishing the 
recursive structure between them. In the proof of this theorem an important role is played by 
the covariant extension (introduced in [6]) of the iV-soliton solution into a new lattice direc- 
tion associated with lattice parameter I = fyv+i- Thus, as in [24], the shift pi — > 'p i denotes 
in what follows the covariant extension of the lattice variables into a direction given by this 
new parameter, thereby exploiting the multidimensional consistency of Q3. Equivalently this 
can be described by the introduction of a new integer lattice variable h associated with the 
parameter k^+i, on which the plane wave-factors (|4,6p depend through the intial values, in the 
usual fashion, i.e. including a new factors with teh lattice parameter 5 replaced by ktv+i and 
the lattice variable n replaced by h. Having this in mind the following theorem describes the 
recursive structure between the elliptic N- and (N + l)-soliton solutions: 

Theorem 6.1. Let be as defined in \5. 3\) and let u^ N+1 ^ be equal to u^ N+1 ^ as defined in 
\5. 3\) , depending on additional parameters rejv+i, cjv+i and additional plane-wave factor pn+i, 
and where all but the latter plane-wave factors pi, i = 1, . . . , N, as well as the discrete exponen- 
tials p(a), p(f3) are replaced by 

^ = e^ + ^- KN+l \ Pl , (i*N + l) , p( K )=e^ M ;- KN+1 \ p( K ), (6.1) 

a(Ki + KN+l ) 0-{K + K N+ l ) 

with k = ±a, ±P respectively. Then u^ N J is related to by the Backlund transformation 

with Backlund parameter £n+i G I\ that is the following equations hold 

Q£ iN+ M N \" (N) M N+1) M N+1) ) = 0, Q$ iN Ju^M N \u^\^ N+1) ) = 0, (6.2) 
in which Qp t N+1) Qq,t N+1 are the quadrilateral expressions of the form given in (|5.1|) 



Proof: The proof of Theorem 16.11 follows a similar pattern as the corresponding proof in 
apart from some refinements which involve the non-autonomous elliptic Cauchy kernel and the 
structure of the elliptic identities involved in the computations. In the first step we break down 
the (N + l)-soliton expression into components associated with the iV-soliton solution. In the 
second step we apply the BT to the A r -soliton solution and in the final step we compare the 
expressions obtained in the two previous steps. 
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. Step # 1. We first establish a recursive structure between the basic objects like S a< p and V a 
between the N- and (JV + l)-soliton solutions. This uses the breakdown of the Cauchy matrix as 
it occurs as the kernel (l + iVf)" 1 in the various objects. First, for the sake of self-containedness, 
we repeat the following lemma: 

Lemma 6.2. The following identity holds for the inverse of a (N + 1) x (N + 1) block-matrix 

(6.3) 



A 




' A-\l + ±bc T A- 1 ) 




c 1 ' 






s 



in which A is an invertible N x N matrix b anc 



are N -component vector column- and 



row-vector respectively , and d is a nonzero scalar, where the scalar quantity s, given by 



d-c T A- v b 



1 


A 


b 


det(A) 


c r 


d 



(6.4) 



is assumed to be nonzero. 

Proof. By direct multiplication, noting that the matrix is invertible if s is nonzero. □ 

Let M (Ar+1) be the (N+l) x (N+l) Cauchy matrix with parameters k\, . . . , kjv+i as defined in 
a similar way as in (|4.4[) . Applying the Lemma to compute the inverse of the matrix 1+M( N+1 > , 
which can be decomposed as above by setting 

A = 1 + MW , b = c N+1 Vt(K N+1 ) X -l +ltK r , 

c T = s t x^ Kn+1 ^^n+i)pn+i , d = l + p n+1 c n+1 ^^(2kn+i) , 

we have 

(i + M^+vy 1 = 



(1 + MM)" 1 + 8- 1 C N+lPN+1 ^(K N+1 )u lt 



N+l U KN+1 



-S 1 C N+1 ^ ! :(k N+1 )u k 

IT— I 



JV+l 



-S i p N+ l^ i (K N+ l) t U KN+1 

where u(-) and u{-) are given in (|4TT71) . Applying the explicit form of 1 + M (Ar+1) ) _1 to the 
definitions (|4.15p to compute the elliptic (N + l)-soliton formulae U^ 1 ^ 1 ^ and Va in terms 
of the iV-soliton ones, we obtain the expressions 



U 



a, /3 



(6.5b) 



v (N + i) = V W + ^ PN+1 c N+1 *l( KN+1 ) (U(%^ 



and subsequently using the definitions (|4.24p we find the following recurrence relations between 
the ./V and N + l elliptic soliton objects: 

(6.6a) 



S 



(N+l) 



C (N) J- C (JV) Q (N) 

a,P ~ ^a,p ~ -CN+lPN+1 ^6c,k n+1 ^k n+1 ,P ' 

W (N + 1) = t ^ 



(6.6b) 



where we exclude values for a and j3 such that a + kat+i = or j3 + ktv+i = (modulo the 
period lattice of the Weierstrass elliptic functions). Furthermore, the quantity s can be directly 
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computed using (j6.4j) with the identifications of the entries given earlier and using the definitions 
(|4.15j) and (14.24[) again, leading to: 

f(N+l) 

s = l + p N+1 c N+ iSjQ} +uKN+1 = j {N) . (6.7) 

Here the latter equality follows from the second equality in ()6.4j) . which expresses s the ratio 
of the r-function f( N+1 ) = det(l + M (Ar+1) ) of the elliptic (N + l)-soliton solution and the 
r-function f( N ^ = det(l + M^ N ') of the embedded elliptic iV-soliton solution. Setting (3 = kjv+i 
in (|6.6ap and using (|6.T|) we subsequently obtain the following identification for the quantity s 

_ >~>a,K N+1 _ k N+1 (( , Q ^ 



da,K N+1 vv k N+1 

which in particular implies that the ratio of S a ^ KN+1 between its N- and (N + l)-elliptic soliton 
value, is independent of the parameter a. 

. Step # 2. Next we apply the BT with the (so far arbitrary) Backlund parameter [ = (/, L) £r a j 
to the elliptic iV-soliton solution defined in (j5.3|) . i.e. we want to solve the system of discrete 
Riccati equations for a new variable v 

Q* l (u^M N \v,v) = , Q* [ (u( N \rt N \v,v) = 0, (6.9) 

relying on the multidimensional consistency of Q3, and the fact that the shifts on the lattice can 
be interpreted as Backlund transformations. To solve the system (|6.9p we reduce the problem by 
identifying two particular solutions, which are obtained by the forward- and backward shifts, in 
the covariant extension, of the known elliptic iV-soliton solutions u^ N \ This implies we need to 
extend all the plane-wave factors by including a lattice direction associated with the parameter 
[, i.e. we set 

po,o(^)=fe" 2 ^)^Vpo,o,o =* m=e 2 ^ J -^ , p(«) = e" 2 ^A- Pi {k) , 



(6.10a) 

in all relevant definitions, where £ K = ^\(k) (A being the uniformising variable associated with 
the point [ on the curve r a ,6 and on the Weierstrass curve). Similarly we have 

p. =e 2v^L^ pi ; £ . =e -2^i^ p .. (6 . 10b) 

Here the elementary discrete shift in the discrete variable h is indicated by a — , in terms of which 
the elliptic iV-soliton solutions depending on three variables, u^^L,, we have 

jriN) (TV) (TV) _JN) 

a n,m,h a n,m,h+l ' —n,m,h 

It is an immediate consequence of the construction behind Theorem 15. 1\ that the following 
equations are satisfied: 

Q^(u^M N \u^M N) ) = , Q*(uW,uWM N \$ N) ) = , (6.11a) 
Q p >W,^), M W,lW) = , Q J[(it (JV) , u (iV) , « (JV) , « (JV) ) = , (6.11b) 
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where (|6.11b|) holds because of the symmetry of the quadrilaterals. Having established the two 
solutions u( N ^ and of (16. 9p we can now find the general solution of that system in the 
interpolating form: 

v= t+T ■ ( " 2) 

in which t is some function to be determined from the following coupled system of homogeneous 
linear equations: 



(6.13a) 



The compatibility of this system is equivalent to the compatibility of the BT (|6.9p as a coupled 
system of discrete Riccati equations, which in turn is a consequence of the multidimensional 
consistency of the lattice equation. 

Using the explicit expressions for the quadrilaterals (|5.ip . eqs. (|6.13p reduce to 

t PuW - {p 2 - l 2 )u^ - L M W w-w 



t PuW - {p 2 - l 2 )uW - LuW w-w 
t QuW - (q 2 - l 2 )uW _ L u(N) m _Q 



(6.14a) 
(6.14b) 



* QuW - {q 2 - l 2 )u( N ) - LuW w-w 

where in the last step we have made use of eqs. (|5.9p . replacing e, respectively 5, by A. Finally, 
using eqs. (|4.34j) with the same replacements, leading to 

_ „ Pa e~^ a W a - £ a e~^ a W a „ p. a e^ a W a + tae-^Wa 

w — w = ~ , w — w = — , (6.15) 

which by setting subsequently a = — A => £ a = yields the expressions 

* = e -(m+vs» ( M W.xW^x t = e - {r , e+fje )x ( qx \ W. X W. X 
t \p„xj W.xW.x ' t \Q-\J W.xW^x ' 

These relations can be simultaneously integrated, yielding the following expression for the func- 
tion t: 

e 2 C (5)A^A_\ / 2 C (e)A^_\ W _ X W_ X . (6.16) 

Using the relation 

W a W p = e a e-^ a S a> p - t-pe^S a>l3 => W X W X = £xe~^ X S x ,x , 

obtained by replacing 5 by A and ~ by — in (|4.25p . and using also (|4.27p with the same replace- 
ments, we can obtain from (16.161) the following result: 

1 a 2 (\)a{2\) p(X) 



t t Po,o(A) 



Sx,x ■ (6.17) 
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Step # 3. We now compare the expressions for the solution obtained from the recurrence of step 
# 1 and the one from the BT of step # 2. On the one hand, the solution from the BT (16.121) 
can be recast into the following form 

v = uW + J_( ¥ W_„W) ^ = tt W_JL L W _W)U (N) , (6.18) 

± -\- 1 L -\- 1 L 

where we have used the relation w — w = L(u^ — u^)/U N ^ , which follows from (|5.9bp . 
setting 6 = e = A. Setting now in (|6.17p the Backlund parameter A = Kjv+i, i.e. I = £n+i = 
(fcjv+ii Kn+i), which implies that the "-shift from now on is the lattice shift associated with the 
lattice parameter fyv+i> and we choose the integration constant to, such that we identify (|6.17p 
as follows: 

l t = c N+ i PN+ iSW +uKN+1 =► « = 1 + |, (6.19) 
comparing the result with (16. 7p . This identification requires the choice 

? = 0-(2KAr + i)cr 2 (KAr + i)/cAr + i/)o.o(KAf+l) . (6.20) 

On the other hand, inserting the recursion relation (|6.6a|) into the expression for the elliptic 
(N + l)-soliton solution we find 

u (N + i) = V /2 («)p 1/2 (/3)5S +1) + --- 

= v /2 («) P i/2 w) to - ^^su ) + • • • 



(AO 1 
S 



(AO 1 

IT ; CN+lPN+1 

S 



A P y\a) P y^)s^ +l s^ +itP+ .. 

—(AT) (JV) (JV) , 2 



where the dots stand for the remaining terms with coefficients B, C, D and with (a, (3) replaced 
by (a, — f3), (—a,/?), (—a, — /?) respectively, and where in the last step we have made use of the 
relations 

W a W x = £ a e-^ a S a , x , WpW x = lpe-^Sp, x , 

which follows from the earlier relation setting (3 = A, and the same with a and [3 interchanged, 
and subsequently taking A = k,n+i- Taking into account the covariant extension of the plane- 
wave factors p(a), p(f3), applying a relation of the form of (|4.7p . we can conclude from this 
computation that 

u (N + l) = U (N) _ C N+ ip N+1 w 2V (N) 

sK N+ i y KN+1 ' ' v ; 

using the expression (|5.6|) and identifying the factor L with K^ + \ through the relation L = 
{l-M-ptp) 1 / 2 for X = k n+1 . 
Finally, noting that 
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we can compare the expressions for v in (|6,18|) and u^ N+1 ) in (|6.21|) and conclude that they are 
the same provided that 

1 1 



-CN+lPN+lS Kr 



and which holds as a consequence of (|6.19|) . Thus, we have established that with the above choice 
(|6.20p of the integration constant to m (|6.16p the Backlund transformed solution coincides with 
the corresponding elliptic (iV + l)-soliton solution given by the Cauchy-matrix objects. This 
then completes the proof of Theorem 16.11 



Theorem 16. II establishes the precise connection between the structure of the iV-soliton solution 
as given by the Cauchy matrix approach, and the way to generate a soliton hierarchy through 
Backlund transforms. What we conclude is that these two approaches coincide up to a subtle 
identification of the relevant constants in the solution. Since, as was remarked in the Corollary of 
section 3 that the soliton solutions of Q3 really live in an extended four-dimensional lattice, the 
precise identification of those constants is of interest, since they contain possibly the additional 
lattice directions. In fact, in [6l[T0] we established the first soliton type solutions for Q4 through 
the Backlund approach, and it is of interest to see how that approach connects to a (yet unknown) 
representation of multi-soliton solutions in terms of a scheme similar to the one set up in this 
paper for Q3. 

7. Elliptic soliton solutions for degenerate subcases of Q3 



In the same way as in the rational soliton case of [24], one can obtain from the elliptic 
iV-soliton solution for Q3 the corresponding elliptic iV-soliton solutions for all "lower" ABS 
equations through limits on the parameters a, (5 (adjusting, where necessary, the coefficients 
A,B,C,D in an appropriate manner). This follows the coalescence scheme of Figure El The 

Q3l ^fQ2] ^[qT 



H3 



H2 



HI 



Figure 3. Coalescence scheme employed to construct elliptic iV-soliton-type so- 
lutions for the degenerate sub cases of equation Q3. 



upper horizontal sequence in this scheme, involving the degenerations of the Q-equations, are 
obtained from performing careful limits of the type (3 — ► a, whilst the vertical limit from Q- to 
H-equations is obtained from the limits a or (3 — > 0. We observe that this degeneration scheme 
roughly follows the pattern of the limits of the NQC equation that were studied in [30] leading 
to the various degenerations of KdV type lattice equations. 
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7.1. Q-degenerations. We first make the following observation which can be independently 
validated: 

Lemma 7.1. The quantity 

Z = F a + (p(ot) - p(£)) Ua,-a , (7.1) 
where U a - a is defined as before, and where we have introduced 

Tl 777 

F a = Fn,m& «) : = C(0 + 3 (C(« " 5 ) " C(« + *)) + y (C(« " e) " C(« + e)) , (7-2) 
ofreys £/ie equation in the following form: 

(p 2 - a 2 )(z - z){z-I)- {q 2 - a 2 ){z - z){z -2) = \(p'{a)) 2 - -^—^] . (7.3) 

4 \p z — or q* — ar / 

Proof. To prove this lemma, let us consider the quantity z as defined in (|7.ip . with F Q given in 
(|7.2p . From the definition we have, using (|4.26p . 



Z-Z = C($ + S)~ C(0 + \ (C(« - 5) - C(« + ^)) + (p(a) - p(0) ?7 a ,-a - (P(«) - P(0) Ua 



—e^ a W a W- a + \ (C(a + <$) + C(« - *) - 2C(a)) 



Pa 2 p(<5) - p(a) 

where the latter form arises by using (|3.8p . Since z remains unchanged under flipping the sign 
of a, because F a = F- a an d U a ^ is symmetric with respect to the interchange of a and 0, we 
also have 



P-a 2 p(<5) - p(a) ' 

Hence, we have 

V 2 p(<5) - p(a)y V 2p(<5)-p(a)y p a P-a 

where the expression on the right-hand side is invariant under interchanges 5 <-> e , ^ <-> ~ , 
taking into account that p a P-a = p{S) — p(a) = p 2 — a 2 . and that r]$ — rjs = n £ — rj £ . Thus, 
we have that 

2 2 /„ 1 p'(a) \ f~ _ 1 p'(a) 

(p — a ) z — z + 



2p(5)-p(a)yV 2p(5)-p(a) 
(q 2 -a 2 )lz-z + l - ^' {a \ -M l-z-l- P ' (a) 



2 p(e) - p(a) / V 2 p(e) - p(a) 

which after some cancellations leads to 

(p 2 - a 2 )(I- z)(?- z) - (q 2 - a 2 )(z- z)$-z) = \{p'(a)) 2 



4 \p 2 — a 2 q 2 — a 2 

which is eq. (|7.3p . □ 
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The identification with the equation Ql (|2.1ap can be made by identifying the lattice param- 
eters 

o _ 1 p'(a) o _ 1 p'(a) (7 A) 

P 2p 2 -a? ' q 2q 2 -a? ' 1 ] 

making it evident that z defined in (|7.ip provides a solution of Ql in that form with the fixed 
parameter 5 2 = 1. The elliptic soliton solutions for the more general parameter case of Ql are 
obtained as an immediate corollary of Theorem 15.11 by considering the limit (3 — ► a. This leads 
to the following result: 



Theorem 7.2. The following formula 

= Ap(a)S aa + Dp{-a)S. 



ii 



+ 2Bz + C 



in which S± 0} ± a , p(±a) are defined as before, where z is given as in Lemma {7A\ and where the 
coefficients are subject to the relation 



AD-B' 



A 



p'(a) 2 

constitutes the elliptic N -soliton solution of the (Qi)a equation, given in the form 
(p 2 — a 2 )(u — u)(u — u) — (q 2 — a 2 )(u — u)(u — u) + A 



9 9 

p z — a z 







(7.5) 



Hence, (|7.5|) leads to a solution of the equation Ql in the form (|2.1ap using the identifications 
of lattice parameters given by (|7.4|) , for generic value of the fixed parameter. 



Proof. Taking the limit = a + 7, and letting 7 — > in (|5.3p the terms with coefficients A and 
D are regular under this limit, whilst the terms with coefficients B and C possess a singularity 
due to the prefactor ^(±a =F P) in the quantities S± a ^p. Thus we have the limits 



M 2 



p 



1/2 



;±a){ l± 7 Fa+ 



+ 2^ 



r 2 a ± | (p(a + 5) - p(a - $)) ± y (p(a + e) - p(a - e))] + 0( 7 3 ) 



X±a,±/3 — ► 
and introducing 



1 

T- 

7 



1 



i-^7 2 p(e) + o(7 J ) 



\^(±2a) 



1 ± 7 r/ ±2 a + tt7 2 (p(2a) - p(£ ± 2a) + t& 2a ) + 0( 7 3 



1 



±7 (P(«) " p(0) + ^ (p'(±«) + p'(0) + 0(7 3 ) , 
X±a,±a + 7 (p(£ ± 2a) - p(a)) + i (p'(£ ± 2a) =f p'(a)) + 0( 7 3 

T Qj/3 = s T X ~} K (l + M)- 1 ^^ (p(/3) -p(t + P + K))r . 



(7.6) 
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c/ Q ,_ Q - 7 r ai _ a + o(7 2 ) , 

-- + (p(a) - p(0) U a ,-a + 7 

7 



oP(£) - (P(«) - P(0 (»7-a ^a-a + T " -a)) 



=> — ► S tt]a + (2a) [r]2a (1 - Xa,aU a ,a) - (p(£ + 2a) - p(a)) 
Assuming the coefficients B and C to behave as: 

B - 50 + 5x7 + - •• , C - C + C l7 + --- , 
where the singularity can be avoided if we choose Co = -Bo- d 

We now consider the limit to the case of the equation Q2. An extension of Lemma 17. II is the 
following statement: 

Lemma 7.3. The following quantity obtained from the expansion of the Q3 elliptic soliton 
formula, namely 



4 



l r 
2 



n 



■m 



p(£) -F 2 a ±^ (p(a + 5) - p(a - <5)) + y (p(a + e) - p(a - e)) 

(p(a) - p(£)) [(^to + Z 7 a) ,-a + 7±t.,T«] 



(7.7) 



o&eys £/ie following relations 

(p? _ a 2) [ (2 - - f ± ) + (4 - 4)(£- i) 

_ a 2 ) [(* - z)(z' ± - f ± ) + (4 - 4)(J- 2) 



4 \^p 2 — a 2 
-(p'(a)) 2 (z + i+?+i)±i(p'(a)) 3 



g 2 — a 2 



± p'(a)p"(a) 
1 



+ 



p 2 — a 2 g 2 — a 2 



(7.8) 



4 + 



Proof. The proof is by direct expansion of the Q3 solution using the expansions u^ N ' 

z + jz+ in the case that we take the coefficients A = D = C = 0, or the expansion v\ N ^ = 
^ + z — 7-z_ in the case we take the coefficients A = D = B = 0, and using the expansion of the 
proof of Theorem 17.21 to obtain the explicit form of the quantities z± ■ □ 

A combination of the quantities z± of Lemma [7. 31 yields the core of the general elliptic soliton 
solution of the equation Q2, as is expressed in the following statement. 
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Theorem 7.4. The following function 

= (p(a) - p(£)) (t„_ + r_„ A ; 



*0 - g 



+ 



2B n 



2B n 



-F«-2(p(a)-p(0)tf«, 



1 /p'» 3(p"(a)) 2 



p/(a) 2 (p'(a)) 2 



R2 



5, 



-p(a)S 



-— ip(-a)5_ a , 

#0 



(7.9) 



where A\, D±, Bq, B\, C\ are constant^ and with the notations as introduced earlier, i.e. 
(|4.15bp . (|4.24p . (|7.6|) and (|7.2|) . and i/ie identifications (|7.4j) . o6e?/s ffoe equation Q2 in the form 

(H3ED. 

Proof. The result is obtained by systematic expansion with respect to the parameter 7 as intro- 
duced in the proof of Corrollary l7.21 The solution (|5.3p of the Q3 equation under the expansions 
given earlier takes the form 

uW = I n _ 1 + n + 7'«i + 0(7 2 ) , (7.10) 
7 

in which the coefficients U-±, uq and u\ are obtained from the expansions of the various ingredi- 
ents of the elliptic iV-soliton solution as given earlier. Here we still have the freedom to impose 
a dependence of the coefficients A, B, C, D on 7 as we want, and we shall chose these such that 
the coefficients tt_i, no are actually constant. This can be achieved by chosing 

B + 7-B1 + 7 2 B 2 , C = -B + 7C1 + 7 2 C 2 



A 



jA\ , D = 7D1 , B = r> -t- 7131 -1- 7--D2 , ^ = —-Do i~ 7L--1 -1- 7-L/2 

in which Ai, Di , Bq, B\, B2, C\, C 2 are all assumed to be of order 1 in powers of 7, and which 
will be specified further lateron. Thus, we obtain u_i = —2B , Uq = C\ — B\ and u\ given by 

ui = (C 2 - B 2 ) + (Ci + + B (z' + + «L) + Aip(a)5 QjQ + D lP {-a)S- a - a , 

where z is given in (|7.ip and in (j7.7f) . Next we employ the same expansions on the Q3 
quadrilateral (15. ip . where for arbitrary coefficients A, B, C, D, A is given by (|5.5p . Thus, on 
the equation we can compute the expansion by noting that up to fourth in powers of 7 we have: 



(P 2 " a 2 ) 



1 1 jM 

2 p L — a 



1 3 

7 

12 ' 



48 



1 2 

4 7 



p"(a) 



+ I (P'(«)) 2 



— a* 2 (p' - 



2^2 



— a 



3 p'(a)p"(a) 3 

2^2 + ' 



O: 



2 (p 2 - a 2 ) 
, 3 (p») 2 



2(P 2 



(p'H) 3 
4 (p 2 - a 2 ) 3 

P / («)P W («) 
(p 2 



,2^2 



+ 



9(p^a)) 2 p» , 15 (p'(a)) 
2 (p 2 



,2^3 



+ 



(p 2 - a 2 ) 4 



+ 0(7 5 ) 



We note that the solution (|7.9|) of Q2 presented in Theorem 17.51 depends effectively on three coefficients, 
A1/B0, D1/B0 and (B\ + C\)/Bq. We could have replaced these coefficients by single symbols, but for the clarity 
of the proof we have abstained from doing that. 



86 



FRANK W NIJHOFF AND JAMES ATKINSON 



and similarly for Q. Inserting this expansion into the Q3 quadrilateral we obtain up to the 
required order in powers of 7 

Qp t q(u,U,U,u) = 

(p 2 — a 2 ){u — u)(u — u) — (q 2 — a 2 )(u — u)(u — u) — f ^7p'(a) + ^-7 2 p"(a)^ (u — u)(u — u) 



I (7 2 (P'(«)) 2 +7V(a)p» + ^7 4 (P"(«)) 2 + ^VW'W 



UU + UU UU + uu 



p 2 — a 2 Q 2 — a 2 



UU + UU UU + UU 



(p 2 — a 2 ) 2 {q 2 — a 2 ) 2 



uu + uu uu + uu 



+ det(A)(p'(a)) 2 



{p 2 — a 2 ) 3 (q 2 — a 2 ) 3 
1 1 



p 2 — a 2 q 2 — a 2 



1+7 



P(a) 2 



+ 



p 2 — a 2 q 2 — a 2 



+-T 
2 



+§P ff (a) 



( p^(a) 

(p 2 — a 2 ) 2 ' (q 2 — a 2 ) 2 ' (p 2 — a 2 )(q 2 — a 2 ) J ' p'(a) 



+ 



+ 



+ 



2 2 2 2 
£r — cr g z — cr 



Inserting the expansion (|7,10p . we find that terms of of order and 1 in powers of 7 (there are 
no negative powers) cancel, provided that 

(det(A)) = \ u \ > (det(A))i = ^u_iu , 

where (det(.A)), indicates the i th order of the expansion of the determinant, and these hold 
identically true for the choices of U-\ and uq given above. The dominant order is then quadratic 
in 7 and yields the following relation 

Q Pi q(ii,S, u, u) = 

= 1 2 \ (P 2 ~ a 2 )(ui - u\){u-y - Si) - (q 2 - a 2 )(ui - 2i)(ui - ui) 



+ 



1 



1 



p 2 — a 2 q 2 — a 2 



1 



(det(A)) 2 - - ( 2u + + u\ + u\ + Si) 



p"'(a) 3(p») 5 



24 - 1 V P'(«) 2 (p'(a)f 



^ 2 -i(p'(«)) 2 



+ 



(p 2 — a 2 ) 2 (p 2 — a 2 )(q 2 — a 2 ) (q 2 — a 2 ) 2 



+ 0(7 3 )- 



(7.11) 



Thus, if Q P) q(n, u, u, u) = as for the Q3 solution we started out with, then the dominant 
coefficient in the expansion in powers of 7 must vanish as well, leading to the conclusion that 
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v = 2— + 



ui , «§-4(det(yi)) 2 l(p"'(a) 3(p"(«)) 2 



6 V P'(a) 2 (p'(a)) 2 



(7.12) 



the function v obeys the Q2 equation in the form (I2.1bj) with the identifications (|7.4j) . □ 

7.2. H-degenerations. The degenerate limits from the solutions of Q3 to those of the El- 
equations are obtained through the limit (3 — > 0. As an immediate corollary of theorem 15.11 we 
have the following result: 



Theorem 7.5. The following formula 

(A + B(-l) n+m ) p 1/2 {a)W a + (C + D(-l) n+m ) p 1/2 {-a)W. 



V (N) _ jji+m 



(7.13) 



in which W± a and p(±a) are defined as before, constitutes the elliptic N-soliton solution of the 
(Hs)a equation, given in the form 



P a (vv + vv) - Q a (w + vv) + 2(p 2 - q 2 )~^r 



. 



(7.14) 



where P a = \J 'a 2 — p 2 , Q a = \J a 2 — q 2 , and where 

A a = p'(a) det(A) . 

Furthermore, the associated solution 

y(N) = -n+m _ B (_yn+rn^ p V\ a )W a + (C - D{-l) n+m ) p l ' 2 (-a)W- a 



(7.15) 



together with v^' are Miura related to the elliptic N-soliton solution w of HI via the relations: 



w — w 



w — w 



P a v - Q a v P a v - Q a v 



V 



V 



P a v — Q a v P a v - Q a v 



(7.16a) 



(7.16b) 



V V 

Proof. The proof is by simple expansions of the elliptic functions in the solution (|5.3p . namely 
by observing that as (3 — > 0, we have 



P 1/2 {±f3) ~ (it)'* 



1 ± (C(e) - nC(5) - mC(e)) + )f (C(£) - n((5) - m((e)) 2 + 0(/3 3 ) 



choosing an appropriate branch of the square root of p(/3), as well as 



(7.17) 



Xa ,(3 - ^-V a + Pp(t + a) + 0(P 2 ) , ^(a + /3) -> * f (a) (l + n a (3 + 0(/? 2 )) , 



XP,K - j ~ VK + /?p(£ + K) + 0(/? 2 ) 



U a ,f3 ^ " Vq) + /3 2 Z a + 0(/3 3 ) , 
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in which 

Z a = s T VK (l + M)- 1 X ~ 1 K r with r) K = C(£ + K) - C(£) - C(K) . (7.18) 

Hence, 

S a ,f3 = *{(Q + P) [1 - XcpU^p] - W a + (3^ i .(a)(7 ]a -Z a ) + 0((3 2 ) (7.19) 
and, thus, we get the limiting behaviour 

U (N) /ft/ + 0(/3 2 ) (7.20a) 

:= { (A - ^(a) [* c (a)(f7a - Z a ) + (C(0 - <(<*) - mC(e)) W a ] 

+ (C - pV2 ( _ a ) [^(-a)^ _ z_ a ) + (C(0 - n((6) - m((e)) W. a ]} , 

(7.20b) 

as given in (17. 13H . In the quadrilateral (15. ip we have the following limits: 

P -> ^P a + 0(J3) , Q + 0(0) and A = p'(a)p'((3) det(A) A a , 

so that the dominant behaviour of the quadrilateral (|5.1|) reduces to the left-hand side of (|7.14|) . 
Finally, with *^(±/3) -» ±/3 _1 (l - ^p(£)/3 2 + 0(/3 3 )) we get the expansions 



Wp = ^(0)(l-Vp) 



l-/3(a r ii )-/3 2 ((a r »7Jr« ) + ^(0) 



^(i/?)^ - ±i(±ir +m (l±/?u; + ^V + 0(/3 3 )) , 
where w is, once again, the solution (|4.33|) of the equation HI, and where 

w' := (C(0 " " mC(e) - s T <* ) 2 - p(f) - 2(s T VK u°) - (s T u ) 2 . (7.21) 

This implies, in particular, that 

UW _> (/r 1 + + tiwW + 0(/3 2 ) , (7.22) 

from which in turn follow that the relations (|5.9|) through these expansion reduce to the relations 
given in (17.16|) as j3 —> 0. □ 

The quantity v' that is obtained in the higher orders of (|7.20ap obeys also some interesting 
relations, which can be easily obtained from (|5.9p by pushing the expansions one step further. 
Thus, corollary, we obtain 

Corollary 7.6. The quantity v' defined in (17.20bp obeys the following Miura type relations 

P a v — Q a v = (w — w)wv , P a v — Q a v = (w — w)wv , (7.23a) 

P a v — Q a v = (w — w)vw , P a v —Q a v = —(w — w)wv , (7.23b) 

where v is the solution of the equation H3 as given in (I7.13P and where w is the solution (14.331) 
of the equation HI and where the parameters P a , Q a are given as in Theorem \7.5 
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It is noteworthy that the quantity w' , i.e. (I7.21|) . that appears in the higher orders of the 
expansions is, in fact, the object that solves the H2 equation. Thus, we can state the following: 

Theorem 7.7. The quantity w' as defined in (|7.2ip . i.e. 

w ' = w 2 - - 2(s T VK u°) - (s T u ) 2 , 

is a solution of the H2 equation, i.e. it solves the set of relations 



vJ + vJ — p 2 = 2ww , vJ + w' — q 2 = 2ww 



(7.24) 

which is a Miura relation (non-auto-BT) between HI and H2, where w is a solution of HI, and 
where w' is a solution of 



(w' + w' - p 2 ){w' + w - p 2 ) = (V + w' - q 2 ){w' + w - q 2 ) 



(7.25) 



which coincides with (|2.4b|) . 



Proof. To prove Theorem 17.71 we proceed in s similar fashion as in section 4. From (|4.18p . using 
To]), we have using the definitions: 

p K e~ VsK (l + M)u° = p K e~^ K r 

= c -^K c -2C(S)K p _ K C 2\C(£)-C(£)]K r 



e VsK ^>s(-K) 



(VK ~Vs)r= (1 + M)p K e~ r)sK - rs T 



o 



u 



where the last step follows by applying (|4.10b|) . Multiplying both sides from the left by (1 + 
M)~ l we obtain the relation 

PKe- VsK u° = (1 + M)~ l n K r- (j] S - (s T « )) u° . (7.26a) 

In a similar way, using (I4.10aj) . we obtain the relation 

P- K e mK u Q = - {l + M)- l ^ K r- (rj s + {s T u )) u° , (7.26b) 

taking into account that from the definitions we have rj_s = —rjg . Multiplying (|7.26ap from the 
left by s T we now compute 



3 T pKe~ IlsK u° = 7? 



^l^ s (K)e^ K u° 



C(0 + C(K) + ((6) - C(e + K)\ u° = -s\ m + r} K )u° 
s T (l + M^tjk r-( m - (S T u )) s T u° , 



and hence we obtain 

- V5 (s T u°) - s T 7} K u° = V VK r- (pg - (s T £t )) (s T u°) . 



(7.27) 
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The same equation is obtained by multiplying (|7.26b|) from the left by s T and performing a 
similar calculation. Using the abbreviations wq := (s T u°) , w\ := (s T r^w ) we can rewrite 
(|7.27|) more conveniently in the form: 

Wl + Wl= TjsWo - 7]SW - W Q W =>• 2(w 1 + Wi) + Wq + Wq + T]g = (w - W + m) 2 ■ 

Using the addition formula for the Weierstrass elliptic function 

V8 = (C(£ + *) - C(0 - CW) 2 = p(£ + + p(0 + p(S) 

and noting that from the definition (|4.33p we have that w — w = 7/5 + wq — wq , we finally obtain; 

(p(0 + wl + 2m - w 2 ) + (p{l) + wl + 2wi - w 2 ) + p(5) = -2w w (7.28) 

which after identifying the expressions between brackets with the object w' in (|7.21|) leads to 
the first part of the Miura transform (|7.24|) setting p(5) = p 2 . The other part of the Miura 
transform follows by the usual replacements, and subsequently the H2 equation for w' is derived 
by applying the trivial identity (2ww)(2wwy = (2ww){2wwY . □ 

An alternative proof of Theorem 17.71 using the expansion (|7.22[) . can be readily given, but 
we believe it is instructive to provide a direct verification based on the form of the solution. We 
note that the Miura relations (|7.24|) were first given in [B]. 
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